Description of orbit propagation models available in CelestLab
A few orbit propagation models exist in CelestLab. Below is the complete list.
Keplerian model (effect of central force)
This is the most simple model, accounting only for the central attraction of the main body.
The main asset of this model is that it can be used for elliptical and hyperbolic orbits.
Secular J2 model (secular effects due to J2)
This model only includes the secular effects on argument of perigee, right ascension of ascending node and mean anomaly due to J2. It is valid for elliptical (including circular) orbits only.
The main asset (in comparison with more accurate models) is that the model can be used for any inclination (including near the critical inclination) and for equatorial orbits.
Eckstein-Hechler model
This model takes into account the effects of the zonal harmonics up to J6 (no tesseral term are considered).
It is adapted to nearly circular orbits (eccentricity less than 0.05 and at most 0.1) and does not work for inclinations close to the critical inclinations (~63.43 deg and ~116.57 deg).
The mean elements include secular and long-period effects.
Note: The CelestLab version of the model includes a small adjustment to the mean eccentricity vector so that the eccentricity of a frozen orbit remains perfectly constant.
References
1) A reliable derivation of the perturbations due to any zonal and tesseral harmonics of the geopotential for nearly-circular satellite orbits, M. C. Eckstein, F. Hechler, ESRO SR-13-1970, ESA, Darmstadt, Germany, 1970.
2) MSLIB FORTRAN 90, CNES, Volume E (me_eck_hech).
Lyddane model
This model, based on Brouwer theory modified by Lyddane, takes into account the effects of the zonal harmonics up to J5.
It is useable for orbits with an eccentricity smaller than 0.9, and does not work for inclinations close to the critical inclinations (~63.43 deg and ~116.57 deg).
There are 2 models derived from Lyddane in CelestLab:
- The original Lyddane model for which the mean elements are orbital elements that vary secularly in time.
- A second Lyddane model for which the mean elements include secular and long-period effects (the mean elements are then comparable with those defined by the Eckstein-Hechler model).
References
1) Solution of the problem of artificial satellite theory without drag, D. Brouwer, 1959, The Astronomical Journal - 64-1.
2) Small eccentricities or inclinations in the Brouwer theory of artificial satellite motion, R.H. Lyddane, 1963, The Astronomical Journal - 68-8.
3) MSLIB FORTRAN 90, CNES, Volume E (me_lyddane).
Note: All the propagation models described above are accessible through the functions: CL_ex_propagate, CL_ex_osc2mean and other dedicated ones.
Other orbit propagation models (available through CelestLabX) are listed below.
STELA propagation model
This model, based on a semi-analytic extrapolation method, takes into account the effects of Earth, Sun and Moon gravity, atmospheric drag and solar radiation pressure.
The effects of the forces are averaged over one orbit (either analytically or numerically) before being integrated. The integration is therefore fast and is suited for long-term propagation.
References
1) Long term orbit propagation techniques developed in the frame of the French Space Act, H. Fraysse et al., 2011, 22nd ISSFD.
The version available in CelestLab is an interface to the original Java code. For more information, see the Stela page.
SGP4/SDP4 model (TLE propagation)
The model (named SGP4 although it is originally refered to as SGP4 or SDP4 depending on the value of the orbit's period) is used to propagate NORAD Two Line Elements (TLE).
The model is based on Brouwer theory.
The effects of atmospheric drag are taken into account through a power density function.
References
1) Spacetrack report #3: Models for propagation of NORAD Element Sets, F.R. Hoots, R.L. Roehrich, 1980.
2) Revisiting Spacetrack Report #3: Rev 2, D.A. Vallado, P. Crawford, R. Hujsack and T.S. Kelso, 2006, AIAA-2006-6753-Rev2.
The version available in CelestLab is an interface to Vallado's C++ code. For more information, see the TLE page.